Nucleation in a sheared Ising model: effects of external field

Simulations using the Forward Flux Sampling method have shown a nonmonotonic de- pendence of the homogeneous nucleation rate on the shear rate for a sheared two dimensional Ising model [R. J. Allen et al, arXiv cond-mat/0805.3029]. For quasi-equilibrium systems (i.e. in the absence of shear), Classical Nucleation Theory (CNT) predicts the dependence of the critical cluster size and the nucleation rate on the external magnetic field. We investigate the behaviour of the sheared Ising model as a function of the external field. At low exter- nal field strength, the same nonmonotonic behaviour holds and the peak in the nucleation rate is remarkably insensitive to the field strength. This suggests that the same external field-dependence holds for the enhancement of nucleation by shear at low shear rates and the suppression of shear at high shear rates. At high field strength, the nucleation behaviour is qualitatively different. We also analyse the size and shape of the largest cluster in the transition state configurations, as a function of the external field. In the sheared system, the transition state cluster becomes larger and more elongated as the field strength decreases. We compare our results for the sheared system to the predictions of the CNT for the quasi- equilibrium case, and find that the CNT cannot easily be used to describe nucleation in the system under shear

Lineage dynamics in growing biofilms: Spatial patterns of standing vs. de novo diversity

Microbial biofilms show high phenotypic and genetic diversity, yet the mechanisms underlying diversity generation and maintenance remain unclear. Here, we investigate how spatial patterns of growth activity within a biofilm lead to spatial patterns of genetic diversity. Using individual-based computer simulations, we show that the active layer of growing cells at the biofilm interface controls the distribution of lineages within the biofilm, and therefore the patterns of standing and de novo diversity. Comparing biofilms of equal size, those with a thick active layer retain more standing diversity, while de novo diversity is more evenly distributed within the biofilm. In contrast, equal-sized biofilms with a thin active layer retain less standing diversity, and their de novo diversity is concentrated at the top of the biofilm, and in fewer lineages. In the context of antimicrobial resistance, biofilms with a thin active layer may be more prone to generate lineages with multiple resistance mutations, and to seed new resistant biofilms via sloughing of resistant cells from the upper layers. Our study reveals fundamental “baseline” mechanisms underlying the patterning of diversity within biofilms

A continuum theory of phase separation kinetics for active Brownian particles

Active Brownian particles (ABPs), when subject to purely repulsive interactions, are known to undergo activity-induced phase separation broadly resembling an equilibrium (attraction-induced) gas-liquid coexistence. Here we present an accurate continuum theory for the dynamics of phase-separating ABPs, derived by direct coarse-graining, capturing leading-order density gradient terms alongside an effective bulk free energy. Such gradient terms do not obey detailed balance; yet we find coarsening dynamics closely resembling that of equilibrium phase separation. Our continuum theory is numerically compared to large-scale direct simulations of ABPs and accurately accounts for domain growth kinetics, domain topologies and coexistence densities

Parasites on parasites:Coupled fluctuations in stacked contact processes

We present a model for host-parasite dynamics which incorporates both vertical and horizontal transmission as well as spatial structure. Our model consists of stacked contact processes (CP), where the dynamics of the host is a simple CP on a lattice while the dynamics of the parasite is a secondary CP which sits on top of the host-occupied sites. In the simplest case, where infection does not incur any cost, we uncover a novel effect: a non-monotonic dependence of parasite prevalence on host turnover. Inspired by natural examples of hyperparasitism, we extend our model to multiple levels of parasites and identify a transition between the maintenance of a finite and infinite number of levels, which we conjecture is connected to a roughening transition in models of surface growth

Kinetic models of ion transport through a nanopore

Kinetic equations for the stationary state distribution function of ions moving through narrow pores are solved for a number of one-dimensional models of single ion transport. Ions move through pores of length $L$, under the action of a constant external field and of a concentration gradient. The interaction of single ions with the confining pore surface and with water molecules inside the pore are modelled by a Fokker-Planck term in the kinetic equation, or by uncorrelated collisions with thermalizing centres distributed along the pore. The temporary binding of ions to polar residues lining the pore is modelled by stopping traps or energy barriers. Analytic expressions for the stationary ion current through the pore are derived for several versions of the model, as functions of key physical parameters. In all cases, saturation of the current at high fields is predicted. Such simple models, for which results are analytic, may prove useful in the study of the current/voltage relations of ion channels through membranes

Computing the local pressure in molecular dynamics simulations

Computer simulations of inhomogeneous soft matter systems often require accurate methods for computing the local pressure. We present a simple derivation, based on the virial relation, of two equivalent expressions for the local (atomistic) pressure in a molecular dynamics simulation. One of these expressions, previously derived by other authors via a different route, involves summation over interactions between particles within the region of interest; the other involves summation over interactions across the boundary of the region of interest. We illustrate our derivation using simulations of a simple osmotic system; both expressions produce accurate results even when the region of interest over which the pressure is measured is very small.Comment: 11 pages, 4 figure